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ABSTRACT 

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet 
powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accom¬ 
panying synchrotron radio emission is produced by the shock interaction between the jet and 
the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, 
however, is the origin of the late-time radio re-brightening, which occurred well after the peak 
of the jetted X-ray emission. Here, we systematically explore several proposed explanations 
for this behaviour by means of multi-dimensional hydrodynamic simulations coupled to a self- 
consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is 
that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one¬ 
dimensional top-hat angular structure. However, a more complex angular structure comprised 
of an ultra-relativistic core (Lorentz factor F ~ 10) surrounded by a slower (F ~ 2) sheath pro¬ 
vides a reasonable fit to the data. Such a geometry could result from the radial structure of the 
super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of 
the ejecta that we infer of ~ few 10^^ erg requires a highly efficient jet launching mechanism. 
Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is 
used to predict the radio light curves for off-axis viewing angles. Implications for the presence 
of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects 
for detecting orphan TDE afterglows with upcoming wide-fleld radio surveys and resolving 
the jet structure with long baseline interferometry, are discussed. 

Key words: radiation mechanisms: non-thermal - galaxies: active galaxies: jets galaxies: 
nuclei - hydrodynamics - radiative transfer 


1 INTRODUCTION 

Radio astronomy is undergoing a revolution in the study of time- 
domain phenomena due to the advent of sensitive wide-held ar¬ 
rays. At met er wavelengths these inc lude LOw Frequency AR- 
ray (LOFAR: Ivan Haarlem et al.ll2013h . the Murchison Widefleld 
Array (MWA; 


Lonsdale et al 


20091) and the Long Wavelength 


Array (LWA; lEllingson et al.lTOO^ At centimetre wavelengths 


(GHz frequencies) a new generation of wide-fleld facilities is 
also being develop ed, such as ApertifAVe sterbork Syn thesis Radio 
Telescope (WSRT: [Oosterloo et al.ll200^ . MeerKAT teooth et al.l 
l2009l) . a nd the Australian Sq uare Kilometer Array Pathfinder 
(ASKAP: Ijohnston et al.ll200^ . However, despite the overall ma¬ 
turity of radio astronomy, surprisingly little is known about what 
astrophysical sources will dominate the transient sky. 

The brightest sources at meter-wavelengths on short 
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timescales are expected to result from co herent processes su ch as 
‘giant pulses’ from Galactic pulsars (e.g. Ijessner et al.ll200^ and 
cyclotron maser emission from brown dwarfs and planets jBergej 
l2002l:lHallinan et al.l2008h . Brighter extra-galactic analogs of such 
events have not yet bee n detected (although see iThomton et al.l 
I2OI3I : [hitler et aHl2014l) . but may accompany rare violent events, 
such as giant magne tar flares 1 L^ubar^ 201^ and the merger of 
binary neutron stars jHansen & Lvutik^ 200 ih . 


At centimeter wavelengths and on longer timescales, in¬ 
coherent synchrotron sources may instead dominate the tran¬ 
sient sky. These generally result from the shock interaction be¬ 
tween rapidly expanding matter from an energetic explosion and 
dense ambient gas . Known examples include radio supernovae 
dWeiler e t al.l l200 2h. off-axis (‘orpha n’) afterglows of gamma- 
ray bursts i Totani & Panaites^l2002h . and other relativistic out- 
flows from the core c ollapse of massive stars (e.g., SN 2009bb; 
ISoderberg et a D I2OIOI) . The transient Swift J164449.3+573451 
(hereafter Swift J1644+57) is the prototype for a new type of rel- 
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ativistic transi ent that will also contribute significantly to the vari¬ 
able radio sky dOiannios & Metzgeill201 ih . 


Swift J1644-H57 wa s characterised by powerful non- 
thermal X-ray emission jBloon^ aUJgO 1 1| ; iLevan et"^ 1201 ll : 
iBurrows et alJIZOlU : IZauderer et al. 2011 1. The long duration of 
Swift J1644-I-57 and a position coincident with the nucleus of a 
previously quiescent galaxy lead to the conclusion that it was pow¬ 
ered by rapid accretion onto the central supermassive black hole 
(SMBH) follow ing the tidal disruption of a star (a ‘ti dal disruption 
event’, or TDE: l^rter & LuminellfToS^ i lReesl l 1^9881), a lthough al¬ 
ternat ive explanations have been proposed (e.g. lOuataert & Kaseij 
l2012h . The rapid X-ray variability suggested an origin internal to 
a jet that was relativistically beaming its radiation along our line 
of sigh t, similar to the bl azar geometry of normal active galactic 
nuclei (lBloometal.ll201lh . Swift J1644+57 was also characterised 
by luminous synchrotron radio emis sion, that brightened g radually 
over the course of several months dZauderer et alj|2011 ). Unlike 
the rapidly varying X-ray emission, the radio emission resulted 
from the shock interaction between the TDE jet and the dense 
external gas surro unding the SMBH dOiannios & Metz^l201ll : 
iBloom et alJl20111) , similar to GRB afterglows. A second jetted 
TDE, Swift J20584-05, with sim ilar X-ray and radio properties to 
Swift J1644-I-57 was reported bv ICenko et alj d2012h . 


Jetted TDEs in principle offer a unique opportunity to witness 
the birth of an AGN, thus providing a natural laboratory to study 
the physics of jet production across a wide range of mass feed¬ 
ing rates. Obtaining a better understanding of Swift J1644-I-57 and 
J2058-I-05 would thus have far-reaching consequences for topics 
such as the physics of relativistic jet formation and super Ed¬ 
dington accretion, the conditions (e.g. distribution of accreting or 
outflowing gas) in nominally quiescent galactic nuclei, and pos¬ 
sibly even the astrophysical origin of ultra-high e nergy cosmic 
rays dFarrar & Gruzinovi 120091 : [Farrar & Pirardl201^ . Despite this 
promise, theory has yet to fully exploit the wealth of data avail¬ 
able for Swift J1644-)-57. Accurate, multi-dimensional models for 
TDE jets are needed to quantify how similar events would appear 
to observers off the jet axis. Such information is necessary to con¬ 
strain the presence of off-axis jets in TDEs detect e d via their quasi¬ 
isotro pic thermal emission teower et al.l 12013l : | van Velzen et al.l 
l2013h or to asses s how jetted TDEs will contribute to future r adio 
transient surveys dGiannios & Metzg^l201ll : lFrail et al.ll2012h . 


Perhaps the greatest mystery associated with 
Swift J1644-I-57 is the unexpected flattening or rebrightening 
of the high frequency radio li ght curves on a ti mescale of a 
few months after its initial rise jBerger et al. I I2OI2I: Fig. [B. This 
behavior contradicts expectations of one-dimensional blast wave 
models that assume energy input proportional to the observed 
X-ray emission, and which instead predict that the high frequency 
emission should decay rapidly once the jet energy is transfered to 
the blast jMetzger et al.ll201^ . This behavior suggests that either 
the temporal or angular structure of the jet is more complex than 
assumed in the basic models usually applied to GRB jets (e.g., 
[Berger et al.ll2012l: IWang et al.ll2014i), or that some key physical 
process has been neglected (e.g., Kumar et al.ll2013h . 


In this paper we model the hydrodynamical interaction be¬ 
tween a TDE jet and the circumnuclear medium (CNM) using 
one and two-dimensional relativistic hydrodynamical simulations. 
Assuming that non-thermal electrons are accelerated at the shock 
fronts, we exploit the simulation results to perform a radiative trans¬ 
fer calculation for the synchrotron radio emission. Our goal is to 
provide a comprehensive study of the structure of the jet that gave 



10" 


10 


■E10 


Figure 1. Radio light curves of Swift J1644+5 7 at observing freq uen- 
cies 1.4, 1.8, 4.9, 15.4, 24.4 and 43.6 GHz from iBerger et al.l 1201211 and 
IZauderer et al.l i2013t) . Also shown are the Swift XRT observations in the 
0.3 — 3 keV band (brown crosses, right axis). Note the broad secondary 
maximum in the radio occuning ~ 6 months after trigger, well after the 
X-ray maximum. 


rise to Swift J1644-)-57, in order to identify which (if any) of the 
proposed models best reproduce the radio data. 

This paper is organised as follows. In ^we review the ob¬ 
servations of Swift JI644-I-57 and summarise previous work on its 
interpretation. §[Bdescribes our model for the jet and the technical 
details of scenarios that we consider. In 0we describe our results. 
§ |5] describes applications and implications of our results, includ¬ 
ing the inferred physical structure of TDE jets powered by super- 
Eddington accretion and the prospects for detecting orphan jetted 
TDEs with future radio surveys. In ^we summarise our results. 


2 A JETTED TIDAL DISRUPTION EVENT 
2.1 Observations of Swift J16444-57 

The two distinct components in the SED of Swift JI644-I-57 in- 
dicated different origins for the X-r ay and radio emission 
teloom et alJl20IlL [Burrows et alj|20li]) . Rapid variability of the 
X-rays placed their origin at small radii close to the SMBH, likely 
from a source internal to the jet itself. Constraints on the brightness 
temperature and variability of the radio emission instead place its 
origin at larger radii, where the relativistic o utflow began to inter¬ 
act w ith the CNM surrounding the black hole dCiannios & Metzgen 
l201lh . Figure[T]shows the observed radio and X-ray light curves of 
Swift 11644-1-57. 

After several days of peak activity, the X-ray luminosity de¬ 
creased as a power law oc in time with a ~ 5/3, consistent 
with the predict ed decline in the fall-back rate M of the disrupted 
star ( lReesll988h . Unless coincidental, the apparent fact that M 
implies that both the jet power Lj tracked the instantaneous accre¬ 
tion rate and that the radiative efficiency and the bulk Lorentz factor 
of the outflow remained relatively constant throughout the evenfl. 

A model for the radio emission of Swi ft 11644-1-57 us i ng the 
first three weeks of data was developed by [Metzger et al.l ( 1201 21 
hereafter MGM12), under the assumption that the X-ray luminos¬ 
ity tracked the instantaneous kinetic power of the jet. This sudden 


* Substantial temporal variation of the bulk Lorentz factor would cause 
equally great changes in the beaming degree of the jet emission, making 
hard to understand why x M for ~ 1 year after the TDE. 
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injection of energy into the CNM produces a forward-reverse shock 
structure, which acts to decelerate the ejecta. The evolution is char¬ 
acterised by two stages: (1) an early phase during which the reverse 
shock crosses back through the relativistically expanding ejecta that 
was released over the initial few-day period of peak activity; (2) a 
later phase during which the relativistic blast wave undergoes self¬ 
similar expansion telandford & McKe3ll976h . 

A striking feature of the early-time radio light curve was the 
presence of an achromatic break at t ~ 10 days (Fig. [T]l. MGM12 
showed that this break could result from the transition between 
stage (1) and (2), and that the temporal slopes of the pre- and 
post-break light curve requires that the CNM density ncnm scale 
with radius as oc with p ~ 2. MGM12 also constrained the 
initial (unshocked) Lorentz factor Fj ~ 10 and opening angle 
0j ~ 1 / F, ~ 0.1 to values rem arkably similar to those of blazar 
jets (e.g. IPushkarev et al.|[200^ . The estimated value of the CNM 
density ~ few cm“^ at radii ~ lO'* cm from the SMBH was some¬ 
what lower than th at measured on a s imilar radial scale from SgA* 
in our own galaxy jBerger et al.l2012h . The reverse shock indirectly 
shaped the early light curves of Swift J1644-1-57 through its effect 
on the hydrodynamics, but the observed synchrotron emission was 
probably dominated by electrons accelerated at the forward shock 
(M GM12fl _ 

I Berger et alj ^20121) presented updated radio light curves of 
Swift J1644-(-57, which showed the surprising radio rebrighten¬ 
ing mentioned earlier (Fig. [TJ. This behavior deviated signifi¬ 
cantly from that predicted by MGM12 for a blast-wave evolv¬ 
ing with constant energy, thus demanding that an additional 
sour ce of energy contribute to t he forward shock at late times 
(e.g. lBarniol Duran & Piranl2013h . Although a temporary rebright¬ 
ening can be created by the self-absorption frequency passing 
through the observers frequency band (MGM12), such a transition 
cannot explain the observed light curves. 


2.2 Proposed explanations for the radio rebrightening 

After the one-zone jet model failed to reproduce the late radio ex¬ 
cess observed for Swift J1644-1-57, several authors proposed expla- 
nations for this apparent incre ase in the energy of the forward shock 
jBarniol Duran & Piranl201^ . We briefly review these possibilities 
here. _ 

iBerger et alJ ilOl'jl proposed that the radio rebrightening can 
be explained as the result of slower material catching up and in- 
jecting energy to the de celerating forward shock at late times (e.g., 
IGranot & Kumaii2006r) . In this model the slower material is ejected 
after the initial ultra relativistic jet stage and it contains ~ 20 times 
more energy than the Fj ~ 10 ejecta. The decline of the X-ray 
lightcurve was believed to reflect the decrease in the jet power cor¬ 
responding to the predicted evolution of the rate of mass fallback 
to the SMBH, indicating that the beaming of the jet emission (and 
therefore the jet Lorentz factor) did not vary substantially over that 
period. It is thus not clear in this scenario what process revived the 
source, resulting in the required slower and more powerful second 
ejection event. On the other hand, the energetic, slower ejecta may 
have been ejected simultaneously with the relativistic component, 
e.g., encasing it (see below ). 

IPe Colle et al.l ( 12012h performed 2D hydrodynamic simula¬ 
tions of the interaction between the TDE jet and its surrounding 


^ With the possible exception of the radio emission during the first ~10 
days after trigger; see below. 


environment. The latter was assumed to include both the extended 
hydr ostatic atmosphere which is produced by retu rning stellar de¬ 
bris jLoeb & UlmeJllfl^ ; iGuillochon et alj|2014h on small radial 
scales (< 10“^ pc) and the CNM of the nuclear cluster on larger 
scales. The cocoon c reated by shocked CN M plays an important 
role in the model of IPe Colle et alj l l2012ll in collimating the jet 
and, potentially, impacting the late-time radio emission. The jet is 
only self-collimated by its environment, however, if its kinetic lu¬ 
minosity is less than a critical value which we estimat e to be ~ 10'*'^ 
erg s“' using the formalism of iBromberg et alj j201 ih . i.e. well be¬ 
low the peak (beaming-corrected) luminosity of Swift J1644-(-57, 
> 10‘'5-‘*f- erg s-‘. 

Radio rebrightening could in principle also result from the jet 
encountering an abrupt change in the radial slope of the CNM den¬ 
sity profile. Such a break might be expected because on small radial 
scales the density profile is expected to approach that of a Bondi ac¬ 
cretion flow onto the SMBH (wcnm ~ while on larger scales 

outside the SMBH sphere of influence the density profile should 
flatten. However, past studies of GRB afterglows show that changes 
in the external density profile produce only minor changes in the 
synchrotron lig ht curve behavior (e.g. iMimica & GianniosluOllI ; 
[Gatetalj201_^ an d refer ences therein) and rarely rebrightening. 

Kumar et alj j2013l) emphasise that electrons accelerated at 


the forward shock are subject to Compton cooling by the X-rays 
from the jet. According to their model, the apparent excess in the 
late radio emission is in fact a deficit in the early-time emission, 
resulting from the synchrotron flux being suppressed due to the 
lower effective cooling frequency of X-ray cooled electrons. Al¬ 
though X-ray cooling will certainly play a role in shaping the early 
radio emission, it becomes less efficient as the jet spreads later¬ 
ally, since a smaller fraction of the freshly accelerated electrons 
(only those within t he opening angle of the original jet) are cooled. 

I Kumar et alj (120131) also consider the effect of X-ray cooling on op¬ 
tically thin emission. As we show here, its effect is negligible for 
the radio bands at which the observed emission is self absorbed. 

A final possibility, which we advocate in this work, is that 
the late radio brightening indicates that the jet responsible for 
Swift J1644-)-57 possessed a more complex angular structure than 
the top hat structure adopted in most previous models. A natural 
alternative is that the jet possesses an ultra-relativistic (Tj ~ 10) 
core responsible for the hard X-ray emission, which is surrounded 
by a wider-angle, mildly relativistic (F^ ~ 2) sheath. Core emis¬ 
sion dominates the early-time emission, whil e the more energetic 
sheath dominates the lat e radio emissio n (e.g. Tchekhovskov et al.l 
I2OI4I ; IWang et^l2014l ; see, however, iLiu et alj|2015h . This sce- 
nario allows both components of the jet to be ejected simultane¬ 
ously during the first weeks following the TDE, but the emission is 
initially dominated by that produced at the external shocks driven 
by the fast core. Such a jet structure is physically motivated if the 
slo w ejecta results from an initial phase of precess ion of the fast 
jet iStone & Loeb|[20T^ ; iTchekhovskov et aijf20T^ or simply as 
the result of the separate fast and slow outflow components ex¬ 
pected to accompany super-Edding t on accretion ( Mineshige et al.l 
I 2 OO 5 I ; ICoughlin & Begelmanl 12014 ; |j iang et 'Zl l2014h . A similar 
fast c ore/slow sheath configu ration is motivated by blazar observa¬ 
tions ( IGhisellini et alj|2005h . indicating that such a structure may 
be a generic feature of relativistic jets. 

The rich data set available for Swift J1644-(-57, now spanning 
over three years, makes it possible to construct a much more de¬ 
tailed model for the afterglow, which explores the effects of late en¬ 
ergy injection, complex angular structure, and X-ray cooling of the 
radiating electrons described above. However, most models to date 
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include only one dimensional hydrodynamics. Multi-dimensional 
effects, such as angular structure and jet spreading or pressure con¬ 
finement hy the shocked external medium, become important as the 
jet decelerates to mildly relativistic speeds. Indeed, since by now, 
~ 4 years after the trigger, even the fastest components of the jet 
have decelerated to F < 2 , modeling the jet as part of a spheri¬ 
cal flow is clearly inadequate. Multidimensional models are also 
necessary to properly address the emission observed by off-axis 
viewers, which peak soon after the jet becomes mildly relativistic. 
Below we describe a series of one and two-dimensional hydrody¬ 
namic simulations that explore systematically the above physical 
effects in order to more accurately model the radio afterglow of 
Swift J1644-H57. 


3 MODEL DESCmPTION AND NUMERICAL 
PROCEDURE 

The numerical code MRGENESIS jMimica et al.l |2007| . l2009ll is 
used to perform ID and 2D relativistic hydrodynamic simulations 
in spherical symmetry. MRGENESIS is a high-resolution, shock¬ 
capturing scheme which uses MPI/OpenMP for hybrid paralleliza¬ 
tion and HDF5 libraries for parallel input and output. 

The TDE jet responsible for Swift J1644-I-57 is modeled as a 
relativistic outflow with a temporally constant Lorentz factor Tj and 
a kinetic luminosity Lj that varies in time with a form motivated by 
the observed (isotropic) X-ray light curve. 


= 


Lj 0 if f < tj 


Tj.o 


- 5/3 


( 1 ) 


if t > ti. 


where fj = 5xl0^s~l week is the approximate dura tion of the 
epoch of peak X-ray activity (e.e.. lBurrows et akllTOllh . The peak 
luminosity Lj_o is related to the total isotropic energy of the jet ac¬ 
cording to E[so - fg Ej(t) df Si 5Lj,ofj/2. Motivated by the isotropic 
energy radiated in the soft X-ray band ~ 3 x 10^^ erg l lUevan et al.l 


1 201 ll ; Burrows et alj2011l) , we require that Eiso > 3 x lO^^erg. The 


true (beaming-corrected) energy of the jet is £ = Eisofg', where 
/fb = (1 - cos0j) is the beaming fraction. In our model for fast, 
narrow jets we assume an initial jet Lorentz factor Tj = 10 and jet 
half-opening angle 9j = l/Tj =0.1. We also consider lower values 
of Tj s = 2 for 0 < Sj_s = l/Fj.s = 0.5 in the case of a wider, slower 
jet, or two-zone models with both fast and slow components. 

The inner radius of the numerical grid, R = rj o, corresponds 
to the location where the jet is injected. The initial size of the ra¬ 
dial grid for ID simulations is /■(,„,_o = 3.2 x 10*^ cm, but this is 
periodically enlarged as the jet propagates to larger radii. This en¬ 
largement occurs in such a way as to keep the numerical resolution 
per unit length constant in the radial direction, i.e. new uniform 
zones are added to the grid each time it is extended. For 2D sim¬ 
ulations the angular coordinate 0 has a fixed range [ 0 , ;r/ 2 ], while 
the jet is injected at £ = rj o across the angular range 0 ^ 0j. The 
ratio of the jet pressure to rest mass density is Pj( 0 /Pj(f)c^ = 0 . 01 , 
so that it is sufficiently small for the initial thermal pressure of the 
jet material to be negligible and, hence, the specific jet enthalpy, 
hj = 1 + s-jd + Pj(t)/p^(t)d = 1 (where ej is the specific inter¬ 
nal jet energy). The jet density pj(t) at the inner boundary is then 
related to the jet luminosity and Lorentz factor according to 


Pj(t) = 


Lj(f) 




( 2 ) 


where vj is the jet velocity. The external CNM is also initialised to a 
cold temperature, with a power-law radial density profile oc 
with k ~ \ - 2. Finally, the jet to external medium effective inertia 
density ratio 77 ^ = Pj/ijF^/jpcnmficnmF^^^) is set to be constant in all 
our simulations and roughly equal to 1000 . 


3.1 ID Simulation Details 

A numerical grid with radial resolution Ar = 1.1 x 10'^ cm ex¬ 
tends from the inner boundary R = ryo = lO'® cm. The latter value 
is taken to be sufficiently small so that our simulations adequately 
follow the initial jet-CNM interactions including the reverse shock 
crossing of the jet. The jet evolution is followed for a total duration 
of ~ 6.7 X 10* s . A total of 1600 snapshots of the jet state are saved 
at regular intervals (every 8.3 x 10“* s) during the first 1.3 x 10* s, af¬ 
terward the time interval between the snapshots is progressively in¬ 
creased, resulting in an additional 785 snapshots being stored over 
the remaining 5.4 x 10* s. The c ollection of saved s napshots is 
used as an input for the SPEV code dMimica et al.l200^ in order to 
compute the synchrotron radio light curves. Appendix lAl provides 
additional technical details on the light curve calculation in SPEV. 


3.2 2D simulation details 

2D models are simulated using the same initial and boundary con¬ 
ditions as in the equivalent ID case. The radial and angular reso¬ 
lution are Ar = 4 x 10*“* cm and A0 = ;r/600 rad, respectively. 
We output the snapshots at exactly the same evolutionary times 
as in ID case (§ 13.11 1. It should be remarked that Ar and AS are 
chosen to make the SPEV calculations feasible, since SPEV must 
solve a ’3 4 - 1 ’-dimensional ray-tracing problem, i.e. the main con¬ 
tributor to the computational complexity of this work is SPEV, not 
MRGENESIS. 


3.3 Non-thermal particles and emission 


Energetic electrons are injected behind the forward shock with an 
energy spectrum /(y) oc y^P, as meant to mimic the effects of e.g. 
diffusive shock acceleration. The forward shock is identified us¬ 
ing th e algorithm described in Sec. 3.1. of ICuesta-Martlnez et al.l 
1 I 2 OI 5 I) . Particle acceleration at the RS is ignored, motivated by the 
lack of evidence for a substantial RS emission component from the 
observed radio light curves of Swift 11644+57 (MGM12; although 
the RS may be contributing to the low-frequency emission during 
the first week post trigger-see below). Tractions and eg of the 
post-shock thermal energy are placed into relativistic electrons and 
the magnetic field, respectively. Optionally, we allow that only a 
fraction of available electrons is accelerated at shocks (we adopt 
= 1 except where indicated in Section |4~^ and in Appendix iBll. 
We adopt p = 2.2, = 0.1 and eg = 10“* as fiducial values of the 

microphysi cal parameters, simi lar to those used to fit GRB after¬ 
glows (e.g., iMimica et al.l2O10|) . 

Injected particles lose energy to synchrotron cooling and adi¬ 
abatic losses, and in some cases we also include inverse Compton 
(IC) cooling by X-ray photons (IC cooling occurs in the Thomson 
regime for parameters of interest). The electron Lorentz factors y 
evolve according to 


dy 

df' 


1 dlnp 4 ccTj 


3 df' 


-y- 


3 nifp 


G'g + “rad) 


(3) 


where f, p, u'g and are the time, fluid density, magnetic field 
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and external radiation energy density, respectively, in the comoving 
frame of the fluid. The value of is only nonzero for models 
including IC cooling f$ l4.1.2b . 


4 RESULTS 

Exploring each of the proposed scenarios for the afterglow of 
Swift J1644+57 (' 32.2b across their full parameter space would re¬ 
quire at least tens of numerical simulations, resulting in a non- 
thermal emission computational time that is currently not feasible 
for two dimensional (axisymmetric) simulations. We thus adopt the 
strategy of first performing a large number of exploratory runs un¬ 
der the assumption of spherical symmetry for the blast wave evo¬ 
lution ( 30 30 . Resulting ID light curves are then compared to 
the observations, allowing us to identify the most promising mod¬ 
els, along with the jet/CNM parameters that best reproduce the ob¬ 
servations. These most promising scenarios are then explored in 
great detail by repeating them using a full two dimensional calcu¬ 
lation ( 34.3b . 


4.1 One-component ID models 

We begin by considering light curves produced by one-component, 
spherically symmetric jets with an assumed (fixed) opening angle 
By Note that although the initial jet opening angle Sj does not factor 
into the hydrodynamic evolution in the one dimensional model, its 
value must still be specified because the angular size of the emitting 
region does factor into the emission calculatiorU 


4.1.1 Fast, narrow jet 

We first consider a fast, narrow jet with an isotropic energy Ejso = 
10^^ erg expanding into a stratified external medium with the radial 
density profile Ucnm = ni8(r/10**cm)“*^, where uig = 1.5 cm“^ and 
k is varied between 1, 3/2, and 2. The jet Lorentz factor and half¬ 
opening angle are assumed to be T^ = 10 and Bj = 0.1, respectively. 

Figurel^shows two snapshots of the ID jet evolution, corre¬ 
sponding to an early time in the jet evolution, shortly before the 
RS crosses the ejecta injected at constant luminosity (t ~ 10^ s; top 
panel), and a later time {t ~ 4.3 x 10^ s; bottom panel). At early 
times both the forward and reverse shocks are still close to the head 
of the jet. Behind the dense jet head, which is produced during the 
initial phase of constant jet luminosity, the density decreases to a 
minimum at ~ 1.5 x lO'^ cm. At yet smaller radii, the density then 
again increases oc (eq. ( 21 ) approaching the point of jet injec¬ 
tion. In the more physical 2D case this low-density tail leads to the 
collapse of the jet channel by the cocoon ( 34.3.1b . 

The relativistic hydrodynamical equations are invariant to an 
overall change in th e normalization of the density and length scale 
(Scheck et al.ll2003) . Equation (|3 shows that an increase in the jet 
luminosity Lj (and hence £^^ 0 ) is equivalent to an increase in the 


^ Contrast this with GRB afterglows, in which the shocked jet material has 
Tshlij > 1 prior to the jet break, i.e., the finite opening angle of the jet 
does not affect the observed flux. However, for TDE jets, Tshdj ^ 1, thereby 
suppressing the obseived flux by a factor of ~ (rshd,)^ with respe ct to that 
from an otherwise equivalent spherical outflow iBloom et alj201 1 I) . Models 
for Swift J1644r-57 that assume a spherically symmetric flow may strongly 
overestimate the afterglow emission of the actual beamed event. 




Figure 2. Upper panel: Number density n (full line), thermal pressure P 
(dashed line) and four-velocity Fv/c (grey line) of the fast, narrow jet at 
the time T ~ 9.6 X 10® s. The reverse and forward shocks are visible as 
discontinuities on the right half of the plot. The density peak marks the 
location of the contact discontinuity in both the top and bottom panels. 
Lower panel: Same as upper panel, but for a later stage in the jet evolution 
(T ~ 4.3 X 10 ’s). By this time the reverse shock has crossed the portion of 
the jet shown, while the region near the front of the jet is gradually entering 
a self-similar evolution. 

normalization of the CMN density uig for a fixed value of rjR. Be¬ 
cause changing £iso/«i8 does make a difference in the calculated ra¬ 
dio emission, we exploit this scaling freedom to perform a limited 
scan of the parameter space without having to perform additional 
simulations. 

Eigure[3shows our results for the synthetic radio light curves 
in comparison to the observations of Swift 11644-(-57 (Eig. [T). We 
set «i8 = 5 cm“^, which sets £iso = 3.33 x 10®^ erg. Each panel cor¬ 
responds to a different radio frequency, with different values of the 
assumed CNM density index k shown with different colours. Thick 
and thin lines show light curves with and without synchrotron self- 
absorption taken into account. As expected, the early-time obser¬ 
vations can be fit reasonably well using the r^^^^ profile (MGM12), 
but the late-time emission is clearly underproduced. The radio 
emission is strongly absorbed at <20 GHz during the first 10 - 15 
days, and the peak at each frequency corresponds to the transition 
from optically thick to optically thin regime. Furthermore, at each 
frequency we can see the achromatic break in the optically thin 
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Figure 3. Radio light curves calculated for our fiducial “fast, narrow” jet 
model (’ ^4.1. It . with Fiso = 3.33 X 10^^ erg, nig = 5 cm“^, = 0.2, 

eg = 2 X 10“^. Thick and thin lines show the light curves with and without 
synchrotron self-absorption respectively. Different panels showing different 
observing frequencies, with line styles showing different assumptions about 
the power-law slope of the radial density profile of the CNM, nciim(r) <x 
where k = I (full lines), 3/2 {dashed), 1 {grey). Shown for comparison with 
circles are the radio observations of Swift J1644-I-57 {Fig.[T). 


light curves (thin lines). As discussed in § |2l this break marks the 
jet deceleration time. 

Figure |4] shows the light curves for a case in which nig is in¬ 
creased to 100 cm“^ (and hence £iso to 6.67 x 10^“' erg). Although 
the late radio emission is now reasonably well fit, the early-time 
emission is somewhat overproduced. However, such narrow jets 
suffer from latteral jet speading effects at rather early stages in their 
evolution that are not probed by ID models. In Sect. 4.3, we explore 
a similar, fast/narrow jet model with 2D simulations demonstrating 
that, as a result of latteral spreading, the jet emission is greatly un¬ 
derpredicts observations on timescales longer than ~ Imonth. We 
have not been able to find any 1-component jet model capable of 
explaining the late radio rebrightening. 


4.1.2 X-ray cooling in the fast jet 

I Kumar et al.l ( 12013h propose that the flat radio light curve in 
Swift J1644-(-57 results from an early-time deficit of the syn¬ 
chrotron emission caused by IC cooling by X-ray photons from the 
jet. We test this idea by running jet simulations including radiative 
losses in equation (|3] according to 


u 


rad 


^A(^obs)r 

Anr'^c 


(4) 


where F is the bulk Lorentz factor of the emitting particles and 
is the X-ray luminosity, which scales in the same way as the the jet 
power (eq. (Tj), i.e. according to 

- 5/3 


Lxitohf) = Av,o max 1 


fobs 


(5) 


where again tj = 5 x 10^ s and we consider several values for the 


Figure 4. Same as Figure |3] comparing radio light curves from a jet with 
higher density and energy (Fiso = 6.67 X 10^^ erg; dashed line) that better 
fits the late-time radio data to the Fiso = 3.33x10^^ erg model from Figure|^ 
(full line) that better fits the early data. In both cases we assume an external 
density profile «ciim ‘x 


overall normalization, Lx,o = 0, 5 x lO'^’ergs^', and 10‘**ergs“*. 
The observer time tabs is related to the simulation time t used in 
equation 0 according to 


rcos 6 

^obs ~ ^ ~ 

C 


( 6 ) 


where r and 9 are the radius and the angle with respect to the line 
of sight of the emitting particle at time t. 

Figure [5] shows the radio light curves including IC cooling, 
calculated assuming a powerful jet {Eiso = 6.67 x lO^"* erg) and pa¬ 
rameters {k = 1.5, fe = 0.1, = 0.002) as necessary to match 

the late radio peak. IC cooling does noticeably influence the op¬ 
tically thin synchrotron emission. However, the effect on the ob¬ 
served, optically-thick emission is almost negligible, even for the 
highest assumed value of L^.q. When the emission is strongly ab¬ 
sorbed the observer only sees the newly injected particles close to 
the external shock, while the particles that have cooled are already 
located further behind and their synchrotron emission does not es¬ 
cape. We thus conclude that IC cooling, within a narrow-jet model, 
cannot explain the late radio rebrightening. 


4.1.3 Slow, wide jet 

To contrast with the fast, narrow jet model (Sec. 14.1.1b we also 
consider a slow, wide jet with T^ = 2 and dj = 0.5 (Sec.|4j HB- Fig- 
urel^shows our results for the slow jet light curves {dashed lines), 
calculated for an external density nig = 50 cm“^ and jet energy 
£iso = 3.33 X lO^'^ erg, in comparison to the equivalent fast jet case 
(full lines). The slow jet can explain the late time peak, but it under¬ 
predicts the early-time emission observed from Swift J1644-b57. 
This is because a slow jet coasts for a long time before appre¬ 
ciably slowing down, since it must first sweep up ~ l/Tj of its 
own rest mass. During this extended coasting phase, the emitting 






























Swift J1644+57: powerful jet with fast core and slow sheath 1 



Figure 5. Radio light curves including the effects of X-ray cooling f 34.1.2t . 
calculated for a jet with £iso = 6.67 X 10^^ erg expanding into a.k = 1.5 
medium for different values of the X-ray luminosity Lx,o (eq. O)- While 
X-ray cooling can suppress the synchrotron emission for an optically thin 
medium (thin lines), its effect on the actual emission is negligible because 
the latter is self absorbed (thick lines). 

area increases and so does the luminosity of the source at bands 
where the emission is self-absorbed. The resulting lightcurve rises 
fast before reaching its peak, with the peak time determined by ei¬ 
ther the thick-thin transition or the onset of the jet deceleration, 
by itself inconsistent with the relatively flat light curve observed 
for Swift J 1644-)-57. The comparison between the light curves pro¬ 
duced by a fast-narrow jet and a slow-wide jet (Fig. EJ is suggestive 
of the fact that different parts of a unique jet, moving at different an¬ 
gles and with different Lorentz factors, could account for both the 
early and the late time light curve of Swift J1644-I-57. 

4.2 Two-component ID models 

The one-component jet models considered thus far can reproduce 
either the early or late time light curves well (Fig. |6), but none 
acceptably reproduce the entire emission. We are thus motivated 
to consider a two-component jet, comprised of a fast core (Fj = 
10 ) with a narrow half-opening angle of 0j = 0.1 rad, which is 
surrounded by a slower sheath (F = 2) from 0.1 to 0.5 rad. Both 
components of the jet are assumed to possess an equal isotropic 
energy E\^ = 4 x lO^'' erg. We employ an external medium with 
a Bcnni oc profile and nn = 60 cm“^. As in previous cases we 
assume ee = 0.1, however we allow to vary modestly about its 
fiducial value 10“^ so as to quantitatively minimize the differences 
among observafional dafa and our synthefic models. 

Figure|7]shows the results of our 2-component ID model. The 
model with eg = 10“^ provides a reasonable representation of the 
well-sampled v > 15 GHz light curve, but it underpredicts the flux 
at low frequencies v < 5 GHz and early times f < 10 days. This 
excess at early times in the observational data appears to require 
either a significant reduction of the optical depth at times r < 10 
days, or an additional emitting region not included in the calcula- 



Figure 6. Compaiing radio light curves from a slow-wide jet (Tj = 2, Bj = 
0.5; dashed lines) to the fast-narrow jet (Tj = 10, Bj = 0.1', full lines). 
Thick and thin lines show models including and neglecting self-absorption, 
respectively. Both jets have Eiso = 3.33 X 10^“* erg and propagate into a 
Kcnm ~ external medium with riig = 50 cm“^. We assume = 0.1 
and eg = 0.001. 



Figure 7. Light curves of a two-component ID jet model ( 114.21 including 
a fast core (Fj = 10; Bj = 0.1) surrounded by a slower sheath (Fj = 2; 
Bj = 0.5, each with the same isotropic energy = 4 X 10^“* erg. Different 
line styles show different choices of the assumed post-shock magnetic field 
strength eg = 10“^ (full line) and 2 X 10“^ (dashed line). 
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tions, and quite likely, a combination of both factors. [Mimica et"^ 
1 I 2 OO 9 I) showed that during the early stages of the jet-CNM inter¬ 
action the reverse shock that decelerates the jet material can be 
particularly powerful. Indeed, for the adopted parameters, the rate 
of dissipation of energy at the reverse shock peaks on a timescale 
comparable to the observed jet duration (~ 10 days in the case of 
Swift J1644-I-57). However, for the fiducial parameters adopted in 
the present paper (including those of the CNM), the emission of the 
reverse shock is strongly absorbed below ~ 10 GHz and, hence, the 
reverse shock contribution alone cannot account for the flux deficit 
of our models at early timefl The lack of a dominant reverse shock 
component at early times is potentially supported by the upper lim¬ 
its of < 2- 10% on the 4.8 G Hz linear polarization on timescales of 
weeks after the BAT trigger ( IWiersema et al.ll2012ll . We note, how¬ 
ever, that if we consider a smaller value of the parameter, the 
synchrotron self-absorption at both the forward and reverse shocks 
is substantially reduced, and even counting only the forward shock 
emission, the flux deficit at early times can be ameliorated (see Ap- 
pendixl^for details). 

The model with eg = 0.002, though slightly overpredicting the 
observed emission, nevertheless provides a good candidate to set 
the input for our full 2D simulations described next. It is desirable 
to slightly overpredict the emission with the ID model since in the 
full 2D case the jet will experience lateral spreading. This leads to 
earlier deceleration of the jet in comparison to spherical models, 
resulting in lower flux at bands where the emission is optically thin 
(see next Section). 


4.3 2D simulations 

Motivated by the improvements in accommodating the light curves 
when considering one-dimensional two component models (Fig.|3, 
we perform 2D simulations of two component jet models with 
Fiso = 4 X lO^'* erg. The jet is injected at the inner grid bound- 
ary (r^.o = 2 x 10*® cm) in the angular range [0,0.5] rad. A fast core 
(Ty^ = 10) is injected between 0 and 0.1 rad, while a slower sheath 
(Tj = 2) is injected between 0.1 and 0.5 rad. The external medium 
into which the jet is propagating has a profile and njg = 60 
cm“^. The jet luminosity changes in time according to Eq.[T] 

The radial resolution of the simulations is Ar = 4 x lO*'* cm 
(Sec. l3.2t . but the radial grid is regularly extended as the jet prop¬ 
agates into the external medium. This is achieved by stopping the 
simulation once the jet head nears the outer radial edge of the cur¬ 
rent grid, extending the radial grid and initialising with the external 
medium. The angular resolution is AS = ;r/600 rad and we always 
simulate the full range [0, ;r/2] using 300 points. 

In order to compute the light curves the snapshots of the grid 
state need to be saved periodically (Sec. l3.lt . The time at which the 
snapshots are written is exactly the same as in ID simulations, thus 
enabling us to directly compare ID and 2D models both in terms 
of their hydrodynamical evolution (Sec. l4.3.Tt . as well as in their 
observational signatures (Sec. 14.3.^ . 
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Figure 8. Snapshots of the hydrodynamic evolution of the two-component 
jet. The time in the upper two panels corresponds to the times of the snap¬ 
shots in Fig.[2]and Fig.[^ In each panel the left half shows the (comoving) 
number density rtth of the thermal fluid, while the right half shows the fluid 
four-velocity Yf^ := Fv/c. Since both fast (Fy = 10) and slow (F.^ = 2) 
jets are injected with the same luminosity per solid angle, «th in the (more 
relativistic) fast jet is initially ~ lower than in the slower jet. This 

can be seen in the upper two panels, where the fast jet is seen as a dark- 
coloured channel close to the symmetry axis. Note that the density color 
scale changes from panel to panel. The spherical boundaries on the left side 
of some panels denote the edge of a radial grid which is periodically ex¬ 
tended, see §[ 3 ] 

4.3.1 Hydrodynamic evolution 

Our simulation starts from sufficiently short distance from the cen¬ 
tral engine so that we can follow all the stages of the jet deceleration 
as the jet sweeps through the CNM. Initially both the slow and fast 
jets are overdense with respect to external medium and propagate 
almost ballistically along their initial opening angle (top panel in 

^ Besides the fact that it is technically very challenging to accurately de¬ 
tect the reverse shock in the hydrodynamic simulations, the reverse shock 
emission is strongly suppressed for the parameters at hand, which justifies 
our choice of not including it in the current models. 
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Figure 9. Same as Fig. [2] but in addition showing the values at the axis of 
the two-component 2D jet simulation (thin lines), and a 2D simulation of a 
single-component fast jet (brown lines). For better comparison, results of a 
ID simulation with the same Ar as the 2D one are shown (see Sec. in). 


Fig.[8}. At this stage, both fast and slow jets move radially and are 
slowed down at the reverse shock (located at distance lO'* cm 
at the second panel in Fig. [^. The reverse shock strongly com¬ 
presses and heats the jet material. The shock front can be clearly 
seen in the right panels of Fig. [8] where the unshocked CNM is 
at rest (black area). Throughout this work, we assume that parti¬ 
cle accelerated at the external shock are responsible for the radio 
emission of Swift J1644-I-57. 

There is a fundamental feature of the structured jet dynamics 
that deserves special attention: the jet is remarkably stable through¬ 
out most of the simulated evolution. An analytical study of the sta¬ 
bility of our configuration is rather complex (and beyond the scope 
of this paper), the reason being that the fast jet is underpressured 
with respect to the slow one. This breaks out a common assump¬ 
tion in the study of relativistic jets, namely that the jet is pressure- 
matched with the surrounding matter. As a result of this, the slow / 
fast jet interface develops a Riemann structure consisting on a pair 
of shocks separated by a contact discontinuity. For the conditions 
we study, both shocks move towards the jet axis, though the one 
shocking the fast jet is faster (and weaker) than the one crossing 
the beam of the slow jet. Across the faster shock, the jet Lorentz 
factor slightly decreases from F/ = 10 to F ~ 9, while across the 


slower shock it decreases almost to the slow jet value (F* = 2). 
In order to estimate the shock front velocities, it is important to 
consider that the fast jet is launched to expand radially at an angle 
Sj = 0.1 and, thus, with a velocity component perpendicular to the 
jet axis Vj_ ~ 6/1 - = 0.1c. Balancing the fast and slow jet 

ram pressure perpendicular to the axis we conclude that the shock 
propagation speed towards the axis is very slow, while the effect 
of the initial (thermal) pressure imbalance between both jets gets 
reduced to a relatively thin layer around the fast jet (see, e.g., the 
narrow red layer surrounding the beam of the fast jet -white shades- 
in the right side of the second panel of Fig.j^. 

Having explained the differences with previous jet stability 
studies caused by the pressure imbalance, we now use the re- 
sults of those works to a nalyze fast/slow jet interface properties. 
[Hardee & HughesI (l2003l) showed that pressure matched, cylindri¬ 
cal jets with a fast-spine and a slow-sheath are more stable than jets 
without slow sheaths. In our 2D model the role of the sheath can 
be played by the slow jet. [Hardee & HughesI ( l2003h show that the 
decrease in the growth rate of Kelvin-Helmholtz instabilities is pro¬ 
portional to the velocity shear Av. In our case, Av = Vfj - Vsj - 0.13c 
(vfj and Vsj are the speeds of injection of the fast and of the slow jets, 
respectively). They also find that the presence of an external mov¬ 
ing shear dampens asymmetric modes and raises the growth rate of 
the fundamental pinch mode. In this regard, we note that, initially 
(T < 100 days) the fast jet does not develop any recollimation in¬ 
duced by the pinching of the slow jet. Later (T > 100 days), the 
bow shock and the reverse shock produced by the slow jet pinch the 
beam of the fast jet only very close to the jet head (y ~ 2.8x lO'^ cm; 
top panel of Fig. Hi. This pinching speeds up the beam of the fast 
jet which, as a result, propagates a bit faster than a pure ID “jet” 
with the same properties. This effect is visible in the top panel of 
Fig.[9l where profiles along the axis of the 2D jet model are com¬ 
pared with an equivalent ID jet model. There we observe that the 
forward shock of the ID jet model lags behind that of the 2D jet 
(located at r =: 2.9 x 10‘^ cm). Also in Fig.|3 we display the axial 
profiles of a 2D jet with the same properties as the fast and narrow 
jet, but without a broader slow jet flanking it (“top-hat” 2D jet). 
Without the stabilising influence on the dynamics of the broader, 
slower jet, the top-hat 2D jet clearly lags behind ID jet model. At 
later times (bottom panel of Fig. H, the situation is reversed, and 
the 2D effects (transfer of transversal momentum) makes the two- 
component jet bow shock to propagate at smaller speed. 

At this point the ram pressure of the fast jet is ~ 3 times larger 
than that of the slow jet and the former “breaks-out” of the latter’s 
bow shock, as can be observed from the third panel of Fig. [8] In the 
subsequent evolution, the initially oblate structure of the outermost 
bow shock of the fast jet becomes increasingly prolate. After T ~ 
20 yr (bottom panel of Fig. the longitudinal sizes of the cavity 
blown by the fast jet and the slow jet cavity become comparable. 

Another important difference we find between a top-hat 2D 
jet model and our two-component jet model is the radial width of 
the bow-shock-to-reverse shock region. As time evolves this radial 
width becomes shorter in the top-hat jet, since there is no dynam¬ 
ical collimation exerted by the broad component. For instance, at 
T ~ 500 days (bottom panel of Fig. and second panel of Fig. [^, 
the radial width of the former region is of ~ 5 x 10*® cm for the top- 
hat 2D jet model, while it is about twice as large in our reference 
sheared 2D jet. This difference increases wifh time, particularly af¬ 
ter the constant injection phases finishes. 

At later times the injection luminosity in both jets decreases. 
The tenuous jet finds itself surrounded by a hot cocoon of shocked 
material. The pressure of the surrounding cocoon causes the jet 
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Figure 10. Same as Fig. 0 but showing the light curves produced by ID 
simulations (full lines, corresponding to dashed lines in Fig.|7) and 2D sim¬ 
ulations of a two-component jet (red lines) for = 0.1, eg = 0.002, as well 
as a 2D simulation of one-component fast and narrow jet (grey lines). Both 
ID and 2D simulations assume Ei^o = 4 X 10^^ erg, an external medium 
with r“3/2 pj.Qfl[g and njg = 60 cm“^. The fast and narrow jet {Ef = 10) 
spans the angle from 0 to 0.1 rad, while the slow and wide jet (Fj = 2) 
extends from 0.1 to 0.5 rad. 


channel to collapse (third panel in Fig.[^. At later times very lit¬ 
tle new injection of energy and momentum takes place. The blast 
wave that bounds the structure of the complex cavity will relax into 
a spherical Sedov solution. While we do not follow its evolution 
until this late stage, the resulting emission (of interest here) is prac¬ 
tically spherically symmetric by the end of the simulation (see next 
Section). 


4.3.2 Light curves 

FigurefTOlshows a comparison between light curves computed from 
ID and 2D simulations with the same initial conditions, as well as 
a light curve computed from a 2D simulation of an one-component 
jet (fast and narrow component only). Since the 2D simulations 
are performed with lower resolution in the radial direction in com¬ 
parison to the ID ones, we have also performed ID simulations 
with reduced resolution (to match the 2D radial grid, dotted line). 
Though, from a dynamical point of view the differences in chang¬ 
ing resolution are negligible, the initial speed of propagation of the 
forward (bow) shoc k of the jet is slightly different (~ lO^^c; 


iMimica et al.l 20091) . which slightly affects the flux arrival times 
and the optical thickness of the overall ejecta. One can see that 
the reduced resolution affects the predicted optically thin emission 
at early times (dotted and full thin curves in Fig. llOb . Flowever, the 
emission at these stages is optically thick and the low resolution 
has negligible effect on the emission (see thick lines). We, there¬ 
fore, expect that the resolution of the 2D simulations is sufficient 
to capture the hydrodynamic evolution of the jet well enough to 
compute the resulting emission (see Appendix A for more details 
on other numerical parameters influencing the emission.) 


At low frequencies the emission is absorbed and ID and 2D 
jets have very similar light curves as long as the jet is expanding 
conically. In the 2D two-component simulations the emission flat¬ 
tens earlier as a result of the jet deceleration and expansion along 
the azimuthal direction. In the ID simulations such expansion is 
not allowed and the jet deceleration takes place more slowly. 

At higher frequencies the 2D two-component jet emission 
peaks earlier and at lower luminosity. The ID model overpredicts 
the optically thin emission by factor of ~ 3 - 5, taking into account 
that our ID treatment already takes into account that the opening 
angle of the jet is finite when calculating the emission. Models 
that assume ID, spherically emitting blast overpredict the resulting 
emission by a large factor (~ 100) for timescales ~ months after 
the TDE. 

At early times, when all models radiate in an optically thick 
regime, the differences among ID and 2D models are negligible, 
and all of them fall at roughly the same distance to the observations. 
However, the one-component 2D jet fails badly to account for the 
observed flux after its light curve peak, when the flux decreases 
much faster than either the two component 2D model or than the 
ID single-component jet models shown in Fig. (3] For brevity, we 
do not show here the light curves obtained with a single component 
2D slow-wide jet, but the effects on the light curve are very similar 
to the ones shown in Fig. |6] and discussed in Sect. 14. 1.3] for ID 
models. Considering that a 2D hydrodynamical modeling is more 
realistic than any ID modeling, we find that the orders of magnitude 
flux deficit in single component 2D fast-narrow jet models at late 
times, together with the flux deficit of single component 2D slow¬ 
wide jet modes at early times is an unequivocal proof that a single 
component jet cannot simultaneously explain both early- and late¬ 
time emission, and argues in favor of a two-layer jet model. 


4.3.3 Origin of non-thermal emission 

As remarked in § 13.31 SPEV models non-thermal particle evolu¬ 
tion by injecting them behind relativistic shock fronts. Figure [TT] 
shows the particle distribution and density for the four hydrody¬ 
namic snapshots discussed in Fig.[8l 

Even though there is substantial dissipation at the reverse 
shocks that can contribute to the emission the first ~ 10 days, we 
only model particles accelerated at forward/external shocks. The 
reason is threefold: (1) the emission from the reverse shock is 
strongly dimmed due to the strong synchrotron self-absorption re¬ 
sulting from the adopted microphysical and dynamical parameters 
of our models; (2) we focus on the month to year timescale evolu¬ 
tion where the reverse shock emission is expected to be less promi¬ 
nent; and (3) computing the reverse shock particle evolution and 
emission is numerically very demanding. 

At early times (top panel in Eig. [ID the particles are con¬ 
fined to a relatively small volume behind the shock fronts and con¬ 
tribute to the observed emission. At later times (bottom two panels 
in Fig. ED a significant fraction of particles stream away from the 
shock front since the distance between the forward shock and the 
contact discontinuity separating the shocked jet material from the 
shocked CNM also grows. It should remarked that we stop injec¬ 
tion when of the electrons becomes less than ~ 2. We have 
checked that our results do not depend on the exact value of the 
cutoff. For r > 10 years the blast wave enters the “deep Newtonian” 
regime, where ~ 1 for all freshly inje ct electrons and the parti- 
cle acceleration spectrum is modified (see lSironi & Gianniodboi^ 
for details). 
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Outflows from the super-Eddington accretion flow are also ex¬ 
pected to reduce the mass accretion rate reaching the SMBH hori¬ 
zon at Ria ~ Rghy ?L factor ~ (RoatlRmY’, where f?out is the outer 
edge of the disk and p ~ i. The fact that our modeling demands 
the sheath possess an energy of ~ few 10^^ erg, comparable to the 
entire available rest mass erg of the tidally disrupted star, 

implies that the accreted matter must have been deposited into the 
disk close to the SMBH. This is consistent with the TDE explana¬ 
tion for Swift J1644-t57 because the low angular momentum of the 
stellar debris causes it to circularise interior to twice the tidal radius 
Ri ~ Ri,{M.IMi,f^ ~ 50(M./10®i?o)^’^’f^g, where = GMjd, 
M, is the black hole mass and the numerical estimate is for a solar 
type star with mass = Mp an d radius R^, = Rq. 

iTchekhovskov et ^ ( l2014h alternatively propose that the 
slo w sheath could be created by an initially precessing fast jet 
(cf. lLei et al.l2013h . which then aligns with the SMBH spin as mag¬ 
netic flux builds up near the hole as mass accretes. In this case the 
wide opening angle of the sheath ~ 0.3 results from the large 
initial misalignment between the spin axis of the SM BH and the 
orbital plane of the disrupted star dStone & Loebll2oT3l . 


5.2 Radio follow-up of TDE flares 

Prior to Swift J1644+57, TDE flares were detected exclusively via 
the thermal emission from their transient accretion disks at soft X- 
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Figure 11. Same as Fig. [8] but showing the number density of the non- 
thermal particles instead of normalised fluid velocity. 


5 APPLICATIONS AND IMPLICATIONS 

5.1 Properties of super-Eddington accretion flows 

Swift J1644-I-57 provides a testbed to study relativistic jet forma¬ 
tion associated with super-Eddington accretion, as may occur in 
TDEs as well as other astrophysical systems such as ultra-luminous 
X-ray sources and gamma-ray bursts (GRBs). Our models show 
that the radio afterglow of Swift J1644-b57 is well explained by a 
structured jet, comprised of a high-Tj core surrounding by a mildly 
relativistic (Tj 2) sheath containing a comparable, or greater, ki¬ 
netic energy. 

The geometry that we infer could simply reflect the generic 
structure of jet/accretion disk systems. If powered by the 
Blandford-Znajek process, the fast jet core could contain those 
magnetic field lines threading the SMBH horizon. The slower 
sheath could, by contrast, represent a mildly relativistic outflow 
from the inner accr etion disk. Powerful disk outflo ws are expected 
theoretically (e.g., iBlandford &^Beg elman| 


by radiation-hydrodynamical i Ohsuga et alJ 


1^99 ), as supported 


200.5h and magneto- 


hydrodynamical (I.Tiang et al.l2014ll simulations of super-Eddington 
accretion. 


X-rays from these events does not in itself rule out the presence of 
a relativistic jet because the X-ray beaming fraction ~ 2/9j > 
100 inferred for Swift J1644-t-57 implies that the vast majority of 
TDEs are viewed off the jet axis. The presence of a jet should never¬ 
theless manifest itself as a late radio brightening once the jet slows 
to non- relativistic speeds, allowing the orphan afterglow to become 
visible IGiannios & MetzgenbOl ih . In our favoured model for the 
jet structure of Swift 11644-1-57, which includes an additional wide- 
angle outflow with a lower, mildly-relativistic initial speed Ty ~ 2, 
radio emission will become luminous even at early times for most 
viewing angles. 

Eollow-up searches have been conducted for radio emis¬ 
sion from approximately twenty TDE candidates discovered ear¬ 
lier by their thermal X- ray/optical/UV emission to search for evi - 
dence of an off-axis jet dvan Velzen et alj20lllBwer et al.l2013h . 
With three ex ceptio ns, most notably IC3599 and RXI1420-1-5334 
dBower et all l2013l) . none show detectable radio emission, in 
some cases to deep limits (Table [T]l. This raises the question 
of whether a jet with properties similar to that responsible for 
Swift 11644-1-57 can be ruled out in the majority of TDEs. 

Eigure [T^ compares these radio upper limits (blue triangles) 
and detections (red circles) as a function of time since the flare to 
our calculation of the off-axis light curves of Swift 11644-)-57 at the 
frequency closest to that of each observation. In almost all cases, 
the model predicted flux exceeds the observational upper limits for 
all viewing angles, constraining the energy of the relativistic jet in 
these systems to be typically an order of magnitude less than that 
in Swift J1644-1-57. The radio detection in RX J1420-5334 is too 
bright to be consistent with an off-axis I1644-l-57-like event, while 
that in IC3599 is somewhat too dim. Nevertheless, future observa¬ 
tions of these events over the coming years can test whether they 
are fading in the predicted manner. The early-time radio emission 
from CSS 100217 also appears much too dim to be consistent with 
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Figure 12. Radio light curves of our fiducial two-component jet model 
( ^4.3.2) for dilferent viewing angles compared to radio detections (red cir¬ 
cles) and upper limits blue triangles) for thermal TDE flare candidates (Ta- 
ble[T). Legends give the viewing angle in radians. 1 GHz radio maps for 0 
and 0.8 rad viewing angles are shown in Fig. 1141 

an ofF-axis jet; this supports the interpretation of lPrake et alj ( l201 ih 
that this event was a nuclear Type tin supernova instead of a TDE. 

The conclusion that most TDE flares are not accompanied by a 
jet similar to Swift J1644+57 raises the question of what conditions 
are necessary to produce a relativistic jet. Does jet formation re¬ 
quire special ‘internal’ conditions, such as highly super-Eddington 
accretion or a high black hole spin, or are particular ‘external’ con¬ 
ditions needed to produce detectable emission, such as a low, or 
high gas density near the SMBH. If the CNM density were lower, 
the timescale for jet deceleration and radio rebrightening would in¬ 
crease somewhat, but this delay is unlikely to be sufficiently long to 
hide the radio emission in these events. On the other hand, if the gas 
density were orders of magnitude higher than in Swift J1644+57, 
then the jet could be choked before reaching large radii. In this 
case the reverse shock is sufficiently strong to slow the jet to sub- 
relativistic speeds over just a few weeks, potentially dimming the 
blast wave emission considerably. 

Powering a TDE jet via the traditional Blandford-Znajek pro¬ 
cess may also require the tidal debris to sweep up a large mag¬ 
netic flux from a fossil accretion disk jTchekhovskov et al.ll2014l : 
iKellev et alj2014lt . However, a net larg e-scale field may no t always 
be necessary to produce a powerful jet dParfrev et al.l201^ . Parfrey 
et al. propose that jets can be driven by magnetic fields amplified 
locally within the disk, a mechanism which is particularly effective 
when the disk rotation is retrograde with respect to that of the black 
hole. Within this interpretation, powerful jets like that responsible 



Figure 13. Radio flux observed at 30 GHz {upper panel) and 3 GHz {lower 
panel) as a function of the angular size 61/2 of the half-light radius for a 
source located at redshift z = 0.1 {thick lines) and z = 0.354 {thin lines). 
Curves are shown for different viewing angles 0, 0.3, 0.6, 0.8 and 7t/2 rad. 
The brown solid bracket shows an approximation of the flux and angu¬ 
lar size required for detection and resolving the source, respectively, with 
VLBI. Time marks the point after which the jet is resolved at 30 GHz for a 
viewing angle of 0 rad (black) and 0.8 rad (green). 


Table 1. Radio Observations of TDE Candidates 


Source 

Ol 

t(fl) 

V 

vLv 

Ref. 


(Mpc) 

(years) 

(GHz) 

(10“ crgs->) 



IC 3599 

88 

21.49 

3 

5. 

1 

RX J1420+5334 

2970 

21.49 

3 

3600 

1 

NGC 5905 

52 

21.91 

3 

< 2.0 

1 

RX J1624+7554 

265 

21.67 

3 

< 12 

1 

RX J1242-11I9 

208 

19.89 

3 

< 8 

1 

SDSS J1323+48 

365 

8.61 

3 

<48 

1 

SDSS J1311-01 

750 

8.21 

3 

< 115 

1 

Dl-9 

1700 

8.0 

5 

< 800 

2 

D3-13 

2000 

7.6 

5 

< 960 

2 

TDEl 

645 

5.4 

5 

< 120 

2 

D23H-1 

910 

4.8 

5 

< 200 

2 

TDE2 

1280 

4.3 

5 

< 590 

2 

PTFlOiya 

1130 

1.6 

5 

< 300 

2 

PSl-lOjh 

820 

0.71 

5 

< 300 

2 

NGC 5905 

75 

6.0 

8.6 

<9 

3 

D3-13 

2000 

1.8 

1.4 

< 1000 

4 

TDE2 

1280 

1.1 

8.4 

< 1650 

5 

CSS100217 

700 

0.3 

7.9 

230 

6 

SDSS J1201+30 

700 

1.4 

7.9 

< 1000 

7 


since first X-ray or optical detection. References: fll lBower'et alj201ft t2 lvan Velzen et all 
l2013t Ol lBade et alll99dlKomossj2002t t41 lGezari et all2008tllB^eil201 It (51 lvan Velzen et all 
I2OI It (fi lDrake et all201 it f7) ISaxton et alHoI^ All quoted upper limits are 5cr. 

for Swift J1644+57 are rare because they are produced only fol¬ 
lowing the disruption of stars with orbits located within the plane 
perpendicular to, and moving retrograde with respect to, the black 
hole spin axis. 

5.3 Prospects for future radio surveys 

iFrail et af] ( l2012h estimate that Swift J1644+57-like events could 
dominate the number of GHz transients on the sky, with a rate 
comparable to or exceeding that of GRB orphan afterglows or 
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T = 0.43 years 



Figure 14. 1 GHz radio maps of a two-component jet for 0 rad and 0.8 rad viewing angle (left and right panels, respectively). Each panel intensity is normalised 
separately. Bottom and upper panels correspond to different times after the TDE (see legends). Note that the jet seen at 0.8 rad is strongly absorbed at early 
times (upper right panel), while it is almost transparent at 3.41 years (lower right panel). 


radio supernovae (cf. iMetzger et al.ll2015h . Such events thus rep¬ 
resent prime targets for upcoming wide-held radio surveys such 
as ASKAP, MeerKAT, WSRT/Apertif, and, ultimately, the Square 
Kilometer Array (SKA). 

Our calculated off-axis light curves can be used to estimate the 
rate of off-axis TDEs detected by future surveys. The rate of jet¬ 
ted TDEs is estimated in two ways. First, the fact that one nearby 
event has been detected by Swift over AT =10 years of moni¬ 
toring out to a redshift z = 0.35 (co-moving volume V 11 
Gpc“^) suggests a rate of ~ l/(KAr) ~ 0.01 Gpc“^ yr“*, which 
given a beaming correction ~ 100 for J1644-1-57 (MGM12) implies 
a rate ‘Rjtde ^ 1 Gpc“^ yr“*. Another approach is to take the esti- 
mate d rate of thermal TDE flares of ~ 10^** - 10~^ galax y"* yr“' 
(e.g. Ivan Velzen & Farr^ 12014 IStone & Metzgej |2014|) . which 
given a local galaxy density ~ 10“^ Mpc“^ corresponds to a TDE 
rate of ~ 100 - 1000 Gpc^ vr~‘ . Assuming that £ 10% of TDEs 
produce jets dBower et an i2013l :l van Velzen et al. 1 12013 I; Fig. [nil, 
this corresponds to a local volumetric rate of T^jtde ^ 10-100 
Gpc“^ yr“' for on-axis jetted TDEs, a factor of ~ 10 - 100 higher 
than the empirical estimate above. Since the evidence for powerful 
jets accompanying other TDE candidates are weak, we adopt here 
the more conservative rate of 'R^yde ~ 1 Gpc“^ yr“* inferred from 
observations. 


A survey with flux sensitivity Fun, can detect events to a lim¬ 
iting distance where Fg is the peak flux of 

the transient at a fiducial distance Dg. If the transient lasts a charac¬ 
teristic duration At, then assuming Euclidean geometry the number 
of detectable TDEs across a fraction fty of the sky at any time is 
approximately given by 


t'tobs 


(4trl3)RjTOEAkyAtDo(FolFu^f'^ 


30 


^jTDE 

lGpc-3yr-‘ 


/ /sky \ / At \ / film 

ia25ji^j\0.5mJy 



(7) 


where we have used Fg = 0.2 mjy for Dg = 10^** cm » 3.2 Gpc and 
V = 1 GHz and have taken A? = 3 yr based on the peak flux and 
duration of our off-axis light curves (Fig. I12t . The flux sensitivity 
is normalised to a value Fiin, ~ 0.5 mJy corresponding to the lOcr 
sensitivity of the planned ASKAP VAST-Deep Multi-Field survey, 
which plans to survey a quarter of the night sky (f^y = 0.25) at 
1.3 GHz to an r ms flux limit cr = 0. 05 mJy with an approximately 
yearly cadence dMurphv et al.ll2013l) . We consider that a high de¬ 
tection threshold (~ lOcr) will be necessary to overcome the large 
number of false positive signals resulting from random statistical 
fluctuations across the large number of observational epochs com¬ 
prising the survey dMetzger et alj20f^ . 

Future surveys similar to VAST Deep could thus detect a few 
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tens of off-axis TDE candidates over the course of a few years, de¬ 
pending on the uncertain volumetric rate !/?jTDE and how represen¬ 
tative Swift J1644-I-57 is of the jetted TDE population. The number 
of TDEs detected by the SKA will be substantially higher than the 
upcoming generation of surveys. However, the greater sensitivity 
of the SKA will extend the detection range to high redshifts z > 1, 
making the detection rate also sensitive to the uncertain redshift 
evolution of the TDE rate per galaxy and the supermass ive black 
hole density (e.e. lHopkins et all200'^ : [^acki et al.l2014h . 


5.4 Resolving the jet structure 

Another critical issue is whether the radio emission from 
Swift J1644+57 or TDEs detected by future radios surveys will be 
resolved by very long baseline interferometry (V LBI) observations 
JOiannios & Metzg^bOl ll ; [Berger et ^l2012li . Eigure [13] shows 
our calculation of flux of the radio emission from our 2D models as 
a function of the angular size 61/2 of the radio emission at 3 and 30 
GHz. In the calculations we ignore the counter jet that is likely to 
be present. As a result, we may underestimate the size of the source 
for misaligned jets. The angular size is calculated by first comput¬ 
ing the barycenter of the intensity of a radio map (see Fig. [T4lfor 
examples of 1 GHz radio maps), and then determining the diameter 
of a circle centered on the barycenter containing half of the flux. 

Although the angular size increases monotonically, the flux 
decreases at late times, making it unclear a priori whether the 
source will remain sufficiently bright to be detected once i t be¬ 
comes large enough to resolve. Following iBerger et alj l l2012h . we 
show for comparison an estimate of the minimum flux required to 
detect the source, which we take as the 5 cr flux sensitivity of the 
EVLA at the closest frequencies to those we consider, 5 and 22 
GHz, respectively. We also assume a best-case angular resolution 
of 0.15 mas and 1.5 mas at 30 and 3 GHz, respectively. 

For an on-axis source at the distance of Swift J1644-I-57 (z 
= 0.354; thin lines), resolving the jet appears to be just pos¬ 
sible at ~ 30 GHz on a timescale of ~ 10 years, i.e. within 
the several years. This motivates ongoing VLBI observations of 
Swift J1644-)-57 over the coming years and decade, although the 
source is unlikely to be resolved by more than a factor of ~ 2 in 
radius before its flux decreases below the detection threshold. 

Prospects are better for a closer source at z = 0.1 thin lines) 
for which we predict the source could be resolved at 30 GHz within 
a year of the TDE. Since off-axis events are more common, there 
is a greater chance that one will be observed substantially closer 
than Swift J1644-I-57 in searches by blind radio surveys ( 95.31 1. For 
every ~ 40 off-axis evenfs af z = 0.354, approximafely one should 
occur at z 5 0.1. Given our estimated detection rate (eq. (7)), future 
surveys like ASKAP could over the course of a few years detect 
a handful of events sufficiently close to resolve the radio structure 
(upon waiting a few years after the peak flux). For well-resolved 
sources, the predicted double-lobed radio morphology could be im¬ 
aged, and the proper motion of the emission peak could be deter¬ 
mined, providing additional constraints on the jet structure. It may 
also be possible to distinguish the structured jet model proposed in 
this work from alternative scenarios. 


5.5 Jet contribution to TDE optical emission 

The bottom panels of Fig.[T2]show jet synchrotron emission could 
also contribute to the optical/UV light curves of TDEs. For on-axis 
events, the peak optical luminosity can reach Lopt ~ 10 '*"'ergs“* 


on timescales of weeks. Optical emission in Swift J1 644-H57 was 
obscu red by ~ 8-10 magnitudes of dust extinction jLevan et al.l 
IMIh . making this prediction impossible to test. Unobscured op¬ 
tical emission was detected in J2058-H05 in coinc idence with the 
final epochs of BAT detection JCenko et alj|201 j) at luminosities 
~ 10'*'* - 10'*^ erg s“' for rest-frame frequencies v ~ 1 - 3 x 10'^ Hz. 
However, the hard spectral slope of this emission, Fy oc A with 
yS ~ 2, was consistent with black body emission peaking in the UV 
rather than non-thermal synchrotron emission. Such bright black- 
body emission_characterises other optically- and UV-selected TDEs 
(e.g.. lArcavi et alj2014l and references therein) and is believed to be 
reprocessed UV/soft X -ray emission from the accretion disk (e.g., 
iGuillochon et alj|2014h . although it is not clear whether all T DEs 
will possess this emission component iStone & Metzgei|[2014h . 

For typical off-axis viewing angles ((fobs <: 0.8), we pre¬ 
dict peak optical luminosities from the TDE jet of Lop, ~ lO'*^ - 
10 "'^ ergs^*, similar to those of supernovae and lasting for several 
months. The chief distinguishing features of TDE jet transients are 
( 1 ) flat, non-thermal spectra and ( 2 ) unusually long-durations of 
several months or longer. Given that the beaming-corrected rates of 
jetted TDEs is at least a factor ~ 10"* times lower than the rate of 
Type la supernovae, it is not surprising that no jetted TDEs have 
yet been detected by optical surveys. 


6 CONCLUSIONS 

This paper explores the evolution and structure of relativistic jets 
created following the tidal disruption of stars by massive black 
holes in the nuclei of galaxies. Our work is motivated by the dis¬ 
covery of Swift 116444-57, the prototypical “jetted TDE” that has 
been the target of extensive and ongoing radio follow-up obser¬ 
vations. The radio afterglow of Swift J1644-)-57 contains valuable 
information on the jet acceleration, emission, and structure. It also 
provides a unique probe of the radial density structure of the CNM 
surrounding otherwise quiescent black holes. 

The unexpected radio rebrightening of Swift J1644-1-57 defied 
the expectations of simple jet models usually applied to GRB af¬ 
terglows, motivating a proliferation of proposed explanations in 
the literature. We have performed a comprehensive exploration of 
these possibilities by means of a large number of hydrodynami- 
cal simulations of the jet-CNM interaction in order to find a jet 
model that can best reproduce the multi-frequency radio data of 
Swift J1644+57. 

After experimenting with a large number of one-dimensional 
models (most of which are not shown here for brevity) with varying 
jet speed, opening angle and CNM density and profile, we conclude 
that no homogeneous, “top hat” jet model adequately reproduces 
the entire light curve evolution. Top hat jets can be made to fit either 
the early radio bump or the late-time plateau, but not both of these 
prominent features simultaneously. We instead favor a jet with a 
hybrid angular structure, comprised of a narrow, fast (F ~ 10) core 
and and a wide slow (F ~ 2) sheath, each containing a compara¬ 
ble amount of energy. Our 2D axisymmetric simulations verify that 
such a model works well in explaining the observations. Such a ge¬ 
ometry could result from the radial structure of a super-Eddington 
accretion flow or as the result of jet precession prior to alignment 
of the jet with the black hole spin axis. 

Radio observations extending from days after the gamma-ray 
trigger to years or longer probe both the initial relativistic jet stages 
of jet-CNM interactions as well the late mildly and non-relativistic 
evolution of the blast. Simple analytical treatments that assume 
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self-similar, spherically expanding blast waves cannot be applied 
to the totality of radio data. Models that assume a spherically sym¬ 
metric situation overpredict the resulting emission several months 
after the TDE by one order of magnitude or more. Models that as¬ 
sume conical jet motion but take into account the finite jet opening 
angle into the emission calculation give much more consistent 
results in comparison to the fully 2D simulations (though discrep¬ 
ancies of order of ~ 3 remain). 

Among the models we have considered in this paper, the ones 
which are closer to the observational data of Swift J1644-)-57 pos¬ 
sesses a combined kinetic energy =: 5 x 10^^ erg in mildly- and 
ultra-relativistic ejecta. Although this precise value is somewhat 
degenerate with our assumed values of the shock microphysical pa¬ 
rameters, a large energy > 10^^ erg is required of any model that can 
successfully account for the late rad io bump (see also lBerger et al.l 
I 2 OI 2 I : [Barriiol Duran & Piranll2013h . This is close to the gravita¬ 
tional energy released by the accretion of a solar mass star, im¬ 
plying that the jet launching mechanism in some TDEs must be 
very efficient. We note that our model underpredicts the early-time 
emission at frequencies <10 Ghz, but this could be corrected for 
by adopting somewhat non-standard values for the microphysical 
parameters (Appendix iBt. 

The dynamical evolution of 2D axisymmetric jets flanked by 
a slower sheath is significantly different than the evolution of (sin¬ 
gle component) “top hat” jet models. The former are much more 
stable against Kelvin-Helmhotz instabilities and display a faster jet 
propagation speed during the first months after the TDE. Further¬ 
more, spine-sheath jets, as the one considered here, remain well 
collimated for much longer time than equivalent “top hat” jets. 
Thus, the aspect ratio of two-component jets is more prolate than 
in one-component jets (the former develop rather spherical cavities 
on time-scales of months). If the jet cavity could be resolved using 
VLBI observations, its aspect ratio could provide us clues about the 
jetted nature of the event. 

The opening angle that we infer for the fast jet core 6j ~ 0.1 
results in an X-ray beaming fraction ff ~ 100. The resulting 
beaming-corrected rate of jetted TDEs, IRjtde ^ 1 Gpc“^ yr“', is 
less than one percent of the TDE rate estimated by other means 
(see, e.g.. lStone & Metzgeill2014l) . This finding is consistent with 
the lack of evidence for off-axis jets in TDEs detected by their ther¬ 
mal accretion disk emission (Fig. m. Special conditions appear 
necessary for the creation of a powerful jet following TDEs. 

Off-axis jetted TDEs represent a promising source for upcom¬ 
ing radio surveys. We present self-consistent calculation of the off- 
axis radio light curves of jetted TDEs, based on a 2D model that 
provides a reasonable fit to the light curve of the on-axis event 
Swift J1644+57. Planned surveys such as ASKAP VAST could 
detect tens of jetted TDEs. The angular size of the afterglow of 
Swift J1644-(-57 at ~ 20 - 30 GHz exceeds the estimated VLBI 
resolution limit of 0.2 mas on a timescale of ~ 10 years, at which 
point the flux density is sufficiently high to be detected. Closer, 
off-axis TDEs detected without X-ray triggers by radio surveys are 
well-resolved on timescales of a year or less, allowing for a detailed 
study of the jet structure. 
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Figure Al. Same as Fig. [5] but showing the result of the models computed 
with lower and higher snapshot frequency. The model with ATsnap = B.3 X 
lO'* s (see text for details) corresponds to the n(r) oc model in Fig.[^ 
The differences in the light curves between the various angular resolutions 
are so small that all models seem to lie on top of each other. 
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APPENDIX A: TECHNICAL DETAILS AND TESTS 

In this section we briefly discuss some technical details of the light 
curve calculation in SPEV. As was already discussed in Sections r3.ll 
and 13.21 the results of a numerical hydrodynamics simulation are 
periodically stored so that at the end of the simulation a collection 
of snapshots can be used to produce the synthetic observations. An 
important parameter is the frequency 1 /ATsnap at which the snap¬ 
shots are being stored. We chose ATsnap = 8.34 x 10'* s because it 
provides a good trade-off between the space-time coverage of the 
jet evolution and the technical limitation^ In Fig. lAll we show that 
numerical light curves computed by SPEV for = 8.34 x 10'* 

s are very close to those computed using AT snap = 4.17 x 10'* and 
1.67 X 10= s. 

Next we study the influence of the angular resolution on the 
light curves. For ID models in this paper we use 510 angular zones 
per radian, i.e. Ng = 51 zones in a narrow jet. Fig. IA2I shows that 
that value is adequate since varying Ng by a factor of a few does not 
change our results. 

= We note that the choice of ATs^ap was in large part constrained by the 
fact tha t size of the intermediate files and the computational time depend as 

AT-„2p (^EoH ). 
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Figure A2. Same as Fig. [3] but showing the result of the models computed 
with lower and higher angulai* resolution. The model with AO = 1.96 x 10“^ 
rad con'esponds to the n{r) oc model in Fig.[3] 
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Figure A3. Same as Fig. [5] but showing the result of the models computed 
with different number of bins Nt,- The model with Nh = 64 con'esponds to 
the n(r) oc r“3/2 j^jodel in Fig.[^ 


Finally, we study the influence of the discreti zation in non- 
therrn al electron momentum space (see Sec. 3.2 of iMimica et al.l 
I 2 OO 9 I for more details). The electron Lorentz factors are initially 
distributed logarithmically in Ni, bins, 


yi ~ Tmin 


s(i-\)KNb-i) 
/max \ 

rmm/ 


(Al) 


where yn,;,, and are initial lower and upper cutoffs. Fig. IA3I 
shows the comparison of the light curves for Nt = 32, 64 and 128. 
We note that the size of the intermediate files scales as ~ Ni, and 
chose Nb = 64 as a good compromise between the file size and 
light curve accuracy. 


APPENDIX B: TWO-COMPONENT LIGHT CURVES 
WITH/e =0.1 

In this section we study the case when the number of elec¬ 
trons accelerated to relativistic energies is a fraction of the 
electrons crossing the shock (see equation 14 in Section 3.3 of 
iMimica & Alovll2012l for the definition and discussion of /,). We 
set = 0.1 and compute light curves for our reference ID two- 
component model (Section and 2D two-component model 
(Section [4.3b , and show the results in Fig. ED It turns out that re¬ 
ducing the value of below the standard one reduces the optical 
depth, thus increasing the flux at all times, but especially at early 
times and at low frequencies. In other words, in the optically thin 
regime the change in by one order of magnitude has very little 
impact on the light-curve, while in the optically thick regime the 
flux is roughly inversely proportional to The peak flux and the 
time when it happens does not strongly depend on In the ID 
models, the late time differences between the cases =0.1 and 
= 1 are small. 

Figure IBTI also shows that with a smaller value of the ID 
model tends to overpredict the early time flux at all frequencies ex¬ 



Figure Bl. Same as Fig. |3] but showing the light curves produced by ID 
(full lines) and 2D two-component (red) and one-component (grey lines) jet 
simulations. We use 6 ^ = 0.1, 65 = 0.002, and show the light curves for Ze = 
I (thick lines) and Ze = 0.1 (thin lines). In all simulations £iso = 4 X 10^“* 
erg. 


cept at 1.4 GHz. However, the decrease in optical thickness due to 
lower is not enough for the 2D models to offset the flux deficit in 
our synthetic multidimensional models at f < 20 days. We note that 
this deficit happens both in the 2D two component jet model and 
in the 2D fast-narrow jet model, but the latter still exhibits orders 
of magnitude flux deficit at late times (thus eliminating a single 
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Table Bl. Difference Dy between the observations and our models. 


V 

(GHz) 

= 1 

II 

p 

= 1 

II 

p 

1.4 

0.02 

1.02 

0.09 

0.04 

1.8 

6.73 

10.20 

2.12 

3.93 

4.9 

9.51 

14.36 

1.73 

4.23 

15.4 

40.37 

40.84 

3.54 

11.46 

24.4 

36.89 

38.23 

4.42 

10.63 

43.6 

25.81 

28.68 

1.52 

8.40 


In the second and third column we show the total en'or (Ea. IBlt between the 
ID jet models computed with parameters = 0.1,1. In the fourth and fifth 
columns, the same en'or is computed for the reference 2D two component 
jet model. The light curves corresponding to all these models are plotted in 
Fig.lBT] 

component model as a suitable explanation for the observed phe¬ 
nomenology). This probably means that the emission originating 
at the reverse shock has to be taken into account in those models 
where the optical depth at early times is smaller (i.e., when <K 1). 

In light of the results shown in Fig. im it is difficult to assess 
by eye which value of 4 is optimal to best explain the observations. 
Thus, we have performed a simple error analysis by computing the 
vertical distance from the data to the model for each observational 
frequency v, defined as e,(y) = data,(v) - modefify), where the in¬ 
dex i runs over observational times and model,(y) is obtained by 
interpolation between the two model points computed at times that 
bracket any given observational data. We define the total error at a 
given observational frequency as 

Z)v : = — V log [cosh(e, (y))], (B 1) 

,=i 

where Ny is the number of observational data points at each fre¬ 
quency y. 

As can be observed in Tab. isn in spite of the fact that the case 
= 0.1 clearly improves the representation of the observational 
data for our reference 2D two component model at y = 1.4 GHz, 
there is no quantitative reason to choose this particular value instead 
of the standard = 1. At the lowest frequency (y = 1.4 GHz) the 
total error is about two times smaller for = 0.1 than for = 1, 
but this difference is compensated by the better behavior of the 
model with = 1 at frequencies y > 1.8 GHz. This quantitative 
analysis reinforces our choice of fiducial parameters (particularly 
the value of f^) to show that our models may accomodate the ob¬ 
servations with moderate and sound variations of the microphysical 
parameters. Taken all together, our models may actually the ob¬ 
servations with a suitable combination of a value of 4 in the range 
[0.1,1] and the contribution of the reverse shock at early times if 
the optical depth is reduced. 
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